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Based on Cohen & Glashow's very special relativity [A. G. Cohen and S. L. Glashow, Phys. 
Rev. Lett. 97 (2006) 021601], we propose an anisotropic modification to the Friedmann- Robertson- 
Walker (FRW) line element. An arbitrarily oriented 1-form is introduced and the FRW spacetime 
becomes of the Randers-Finsler type. The 1-form picks out a privileged axis in the universe. Thus, 
the cosmological redshift as well as the Hubble diagram of the type la supernovae (SNe la) becomes 
anisotropic. By directly analyzing the Union2.1 compilation, we obtain the privileged axis pointing 
to (l,b) = (242° ± 44°, -42° ± 23°) (68% C.L.). This privileged axis is close to those obtained by 
comparing the best-fit Hubble diagrams in pairs of hemispheres. 
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1. INTRODUCTION 

The cosmological principle is one of foundations of the 
standard cosmological model, i.e., the ACDM model 
It says that the universe is statistically homogeneous and 
isotropic at large scale. The ACDM model is well con- 
sistent with present cosmological observations, such as 
the Wilkinson Microwave Anisotropy Probe (WMAP) Q 
and the Sloan Digital Sky Survey (SDSS) Q, etc. How- 
ever, there exist challenges for the standard cosmological 
model (see review in Ref. 0] ) , such as the large-scale cos- 
mic flows [5JI , the alignment of low multipoles in the CMB 
spectra [6r|8| , the large-scale alignment of the quasar po- 
larization vectors [§| . One of their resolutions refers to a 
privileged axis at large scale in the universe fioj . 

The type la supernovae (SNe la) have been used to 
search for the possible anisotropy of the universe (Tl - 22 1 . 



since they were employed to discover the cosmic accel- 



eration 



23l |24| . Especially, Antoniou & Perivolaropou- 



los used the hemisphere comparison method to ana- 
lyze the Union2 data [25J and found a direction (I, b) = 
(309 o l3? ,18°Iio°) for the 

maximum accelerating ex- 
pansion of the universe. Similarly, Cai & Tuo (26j 



found a preferred direction (l,b) = (314°I^„, 28°l33o) 
for the cosmological deceleration parameter. Most re- 
cently, Kalus et al. j2?| also found that the highest 
expansion rate of the universe towards the direction 
(1,6) « (-35°, -19°) (95 % C.L.). 

A privileged axis, as is mentioned above, may account 



for the anisotropic phenomenologies of the observational 
astrophysics and cosmology. Actually, the issue of privi- 
leged axis has been studied extensively in the very special 
relativity (VSR) [28|. The Finslerian spacetime struc- 

1-6 t 

ture dr = {-q^dx^dx^) 2 (n a dx' T ) w as p roved to be 
invariant under the DISIM(,(2) group [29J. There is a 
preferred axis n CT in the above line element. It could 
characterize the anisotropy of the flat spacetime, which 
leads to the Lorentz invariance violation (LIV). Locally, 
the Randcrs metric was proved to possess the symmetry 
of the group TE{2) ^E^. The group TE{2) is a semi- 
product of T(4) and E(2) The group E(2) denotes 
a proper subgroup of the Lorentz group with three gen- 
erators [28|. Therefore, Randers spacetime could possess 



local symmetry of the generic VSR. 

In this paper, we propose a Randers line element 
(structure) [33( with local symmetry of the generic VSR 
to estimate the possible anisotropy of the Hubble dia- 
gram of SNe la. A modified cosmological redshift for- 
mula is presented. The redshift is direction dependent 
and refers a privileged axis. This may imply that the 
universe undergoes an anisotropic expansion. Next, we 
show a modified luminosity-distance vs. redshift relation 
for the SNe la. It is certainly anisotropic. The Union2.1 
compilation 34| is used to constrain the direction of the 
privileged axis. The rest of the paper is arranged as fol- 
lows. In section 2, the anisotropic Hubble diagram is 
showed in a Randcrs spacetime. In section 3, we investi- 
gate the Union2.1 dataset of the supernovae to constrain 
the level of the anisotropy of the universe. Conclusions 
and discussions arc listed in section 4. 
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2. THE ANISOTROPIC HUBBLE DIAGRAM 

In the ACDM model, the cosmic spacetime is de- 
scribed by the spatially flat Friedmann-Robertson- 
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Walker (FRW) line element [1] 

3r = ^dt 2 - a 2 (t) [dr 2 + r 2 (d9 2 + sin 2 9 dip 2 )} , (1) 

where a(t) denotes the scale factor of the universe at the 
time t. The cosmological redshift z is given by 



1 + z(t) 



a(t ) = J_ 
a(t) ~ a{t) 



(2) 



where we set a(to) = 1 for today. The redshift de- 
scribes the expansion rate of the universe between t and 
to- The overline denotes physical objects in the FRW- 
Riemannian spacetime. 

The FRW line element should be modified when there 
exists a privileged axis in the universe. We postulate that 
an extra closed 1-form is added into the FRW structure. 
The 1-form singles out a privileged axis in the universe. 
The spacetime structure becomes 



dr = dr 



b fl (x)dx 11 



Jv\\ 



(3) 



0. Here 



^ v denotes the covariant derivative of b(x) with respect 



where the 1-form is closed 
6 

to the FRW metric [35j. Actually, the spacetime struc 
ture fl3J) belongs to Randers type [3j]. This is a class of 
Finsler spacetime (35j. The 1-form could be viewed as 
an arbitrarily oriented electromagnetic 4-potential in the 
universe. It may be the relic of a primordial magnetic 
field at large scales 361-3 



In the FRW-Randers spacetime, the cosmological red- 
shift z could be derived via resolvin g t he Finslerian 
geodesic equations, which are given by [40j 



d 2 x° 
~d^ 



dr 2 



dx l dx 3 dx 



Sijaa— — — 
dr dr 

a dx° dx^ 



25. 



dr dr 



dr 
dx i 
~d^ 



I 



dx 

dx 
d7 



0, (4) 
0, (5) 



where / (x, g~) = b v \\^—^r/F and the dots over a de- 
note j^o. The above geodesic equations have solutions 

as a^- cx Ji and a 2 ^- oc J x , where J x = 1 — b^ [ioj . 
Thus, the redshift z could be written as 



1 



1 



a(t) 



(6) 



where the unit 4-vector p denotes a light-like direction 
towards each SN la. It is anisotropic since there is a 
privileged axis b /1 (x) in the above formula. 

When the privileged axis is a spatial 3-vector, the mod- 
ified redshift z could be rewritten as 



1 



l + z(t,p) = — [l-D(n-p)] 

a(t) 



(7) 



where n is the unit 3-vector for the spatial components 
of 6, and D denotes the magnitude which is smaller than 



one. To simplify our following discussions, we choose the 
spatial 3-vcctor n as the third spatial axis. Thus, the 
redshift z becomes 



1 



1 + z(t,cos6) = - [1 - D cos 



(8) 



where 8 denotes the angle between n and p. This formula 
shows clearly that the cosmological redshift is direction 
dependent and anisotropic. Thus, the universe undergoes 
an anisotropic expansion, which takes a dipole form. 

The Hubble diagram would be anisotropic and have a 
corresponded privileged axis which exists in the Randers 
spacetime. We will discuss this proposition in the fol- 
lowing paragraphs. For the null geodesic, the spacetime 
line element vanishes, i.e., dr — 0. In the FRW-Randers 
spacetime, the Randers structure §3§ should vanishes. 
Thus, we obtain 



dr 



a 2 dr 2 = {-DcosOdrY 



(9) 



where we use the polar coordinates and n is set again as 
the third spatial axis. Here the angular coordinates are 
discarded since they are irrelative to the distance-redshift 
relation. Finally, the above equation could be simplified 

as 



dt 



Va 2 



D 2 cos 2 6dr . 



(10) 



The privileged axis emerges again, which affects the prop- 
agation of the cosmic light. This result is different from 
that of dt = a(t)dr in the ACDM model. The modifica- 
tion is quadrupolar in Eq. (|10p . Obviously, the cosmic 
light would propagate with different speeds in different 
directions in the space. However, the quadrupolar effect 
is a second-order effect. It just affects the propagation of 
the cosmic light slightly. 

In the observational universe, the luminosity distance 
d L of a SN la is given by [4l| 



(1 + z)r 



111) 



where r denotes the comoving distance between us and 
the SN la. We could obtain the comoving distance r by 
integrating the equation (|10| . Therefore, the luminosity- 
distance vs. redshift relation can be obtained as 



,i(z) 

(1 + 2) / B{a,6)dt 



where 

B(a,0) = (a 2 + L> 2 cos 2 i 
The Eq. (|12p could be rewritten as 



-1/2 



1 + z 



B(a,i 



da 



aH^Q, m a 3 + SI/ 



(12) 



(13) 



(14) 
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where fl m = (8irG/3HQ)p m o and f^A = A/3-ffg, p m o is 
the critical mass density, A is the cosmological constant, 
and H is the Hubble constant. In the derivation of (fT4")l . 
we have used the conventional Friedmann equation and 
the definition H = (da/dt)/a in the ACDM model as a 
first-order approximation. 



3. NUMERICAL RESULTS FROM SUPERNOVAE 

We define a Cartesian coordinate system for the unit 3- 
vector p corresponding to each SN la with the equatorial 
coordinates (a, 5), in which p is given as 

p = cos(<5) cos(ev) i + cos((5) sin(a) j + sin(i5) k , (15) 

where i, j and k are basis 3-vectors. Suppose n is given 
by 

n = cos((5o) cos(ao) i+cos(<5o) sin(a ) j+sin(5 ) k , (16) 

then the cosine of the angle 8 between these two vectors 
as 

cos6> = n • p . (17) 
Substituting the relations (JSj) and ([T7f into (fT4l) . one has 
H Q d L 



1 + z 



, dz'i{± + z>) 

B{z ,6)— _ == 

'n. ( t a 7 ) ' + o A 



(18) 



where 



A 



B{z',e) = 



1 - D cos e 



A 



1 + z 1 



n-1/2 



D cos 2 9 



(19) 
(20) 



These relations refer to the anisotropic Hubble diagram 
of SNe la. The direction dependence is obvious in the 
above equations since different modifications are intro- 
duced to the distance-redshift relation of SNe la in dif- 
ferent spatial directions. Once again, we see that the 
1-form in the Randers structure leads to the privileged 
axis in the universe. In the ACDM model, D vanishes. 
Thus, A and B reduce back to 1 and 1 + z, respectively. 

We base our numerical study on the Union2.1 dataset 
[34j| . The Union2.1 compilation consists of 580 SNe la. 
Only 389 of them have available equatorial position data 
in the database of the Central Bureau for Astronomical 
Telegrams (CBAT) [42]. Thus, we perform a least-% 2 fit 



to the data of these SNe la to determine the privileged 
axis (cvrj, So) and the geometrical parameter D (viewed as 
a constant for simplicity). The \ 2 statistic in our fit is 



[fMh(zi, on, Sf, a , S , D) - p hs{zi)f 



i=l 



(j^ziY 



(21) 



where p, t h(zi, cti, Si; otQ, 8q, D) is the distance modulus 
calculated from (fT8|) by using the definition n = 
51og 10 [c?L(Mpc)] + 25. p h s {zi) and a^(zi) respectively 
denote the observational values of the distance modulus 
and measurement errors, which are obtained from the 
Union2.1 compilation. As the first-order approximation, 
the parameters 17a = 0.73, f2 m = 0.27, and the Hubble 
parameter h = 0.72 are used in our fit. In the FRW- 
Randers spacetime, D is around ~ 0.1, which leads to 
only small perturbations to the equation (JTHJ) • 

The best-fit results are given as D = —0.13 ± 0.05, 
and (a ,5 ) = (68° ± 44°, -39° ± 23°) in the equato- 
rial coordinate system or (l Q , b ) = (242° ± 44°, -42° ± 
23°) in the galactic coordinate system, for a minimum 
value of X 2 m i n — 0.99. All the results are presented 
in a sense of 68% confidence level (C.L.). The dis- 
tance modulus vs. redshift relation of the SNe la is 
shown in Fig. [1] This privileged axis obtained by di- 




4i- 



0.6 0.8 1.0 
redshift z 



FIG. 1: The distance modulus vs. redshift relation of the SNe 
la in the FRW-Randers spacetime. The (black) dots and er- 
rorbars denote experimental data which comes from Union2.1 
compilation [33] ■ The (red) triangles denote theoretical pre- 
dictions for the SNe la in our model. 



rect analysis is close to those obtained by comparing 
the best-fit Hubble diagrams in pairs of hemispheres, see 
Fig. [5] For example, Antoniou & Perivolaropoulos flfl ] 
showed (Z , b ) = (309°^?° , 18°± 1( f„), Cai & Tuo [H got 
(l ,b ) = (314°l?g°., 28° 



), and Kalus et al. [27[ ob 



tained (Z ,6 ) w (-35°, -19°) (95 % C.L.). 

4. CONCLUSIONS AND REMARKS 

In this paper, we proposed an anisotropic modifica- 
tion to the FRW line element. The modified line element 
refers to the Randers spacetime, which possesses the local 
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(b=90°) 



(l=60°) 6 




•o"(b=0°) 



Riemannian. Thus, it might be a reasonable candidate to 
account for the privileged axis and the anisotropic prop- 
erties of the universe. 

We thank useful discussions with Jian-Ping Dai, Yun- 
guo Jiang, Danning Li, and Hai-Nan Lin. This work is 
supported by the National Natural Science Fund of China 
under Grant No. 11075166. 



(b=-90 u ) 



FIG. 2: The privileged axis in the galactic coordinate. The 
star denotes the direction of the privileged axis obtained by 
the direct analysis in this paper. For contr ast, we also show 
directions for other privileged axes [lol l26l . [23 obtained via 
the hemisphere comparison method. 



symmetry of Cohen & Glashow's VSR. The local symme- 
try involves the group TE(2). The Euclidean group E(2) 



Jy, -12 



To = K„ 



J x 



contains three generators T\ 
and J z , where Ki and Ji (i = x, y, z) denote the genera- 
tors of boosts and rotations, respectively. The former two 
generators form a two-parameter group of translations in 
the x — y plane. Thus, the local FRW-Randers spacetime 
is cylindrically symmetric. However, the parity violates 
in the z-direction. Otherwise, the E(2) would be enlarged 
to the Lorentz group (28|. Actually, the 1-form b^dx^ in 
the Randers structure changes its sign under the direc- 
tion reversal dx^/dr — > —dx^/dr. This reveals that the 
Randers structure is asymmetric. On the other hand, the 
indicatrix of the FRW-Randers spacetime is quadratic, 
while the center of the indicatrix is displaced in the in- 
direction (33l | . The above two accounts imply a privileged 
axis in the FRW-Randers spacetime. The existence of the 
privileged axis contradicts with the cosmological princi- 
ple and implies a statistically anisotropic universe. We 
extracted the direction (fo, &o) = (242° ±44°, -42° ±23°) 
for the privileged axis b % , based on the Union2.1 com- 
pilation of the SNe la. The direction is close to those 
obtained by comparing the best-fit Hubble diagrams in 
pairs of hemispheres. 

It is noteworthy that the Randers spacetime belongs 
to Finsler geometry [35| . Actually, Finsler geometry gets 
rid of the quadratic restriction on the spacetime struc- 
ture (line element) j43|. It is a natural generalization 
of Riemann geometry and includes Ricmann geometry 
as a special case. Finsler spacetime could depend on 
certain preferred directions of the spacetime background 
29, 44 49]. This could also be revealed via the isomet- 
ric transformation US Ed, GUI. There are no more than 

Killing vectors in the d-dimensional Finsler 
spacetime |50j. Otherwise, Finsler spacetime becomes 
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